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ABSTRACT 


The  propagation  of  cylindrical  dielectric  waveguide 
inodes  near  cut-off  and  far  from  cut-off  are  considered.  The 
relative  amounts  of  E2  and  Hz,  and  the  transverse  components 
of  the  fiel-  are  determined  for  both  sets  of  hybrid  modes. 

With  the  radial  dependence  of  the  z-conponents  of  the  field  in 

the  central  dielectric  given  by  J  (ur/a) ,  the  transverse 

(  “ 

components  far  from  cut-off  are  given  by  Jjj^tur/a),  where  u  is 
a  parameter  found  from  the  boundary  conditions  and  which  fixes 
the  scale  of  the  Bessel  function  relative  to  the  boundary  r=a. 
The  two  \alues  n+1  and  n-1  correspond  to  the  two  sets  of  modes. 
The  designation  of  the  hybrid  modes  are  discussed.  Field  plots 
for  the  lower  order  modes  are  given. 
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LIST  OF  ILLUSTRATIONS 

Construction  to  show  the  equivalence  of  the  waveguide  condition 
that  at  cut-off  Vph  =  c/n2  with  the  geometrical  optics 
conditions  that  propagation  occurs  only  if  the  angle  of 
incidence  of  the  wave  on  the  fiber  wall  exceeds  the  critical 
angle  for  total  internal  reflection.  The  wave  normal  is 
given  by  p.  S.  and  S  are  two  equiphase  surfaces 
separated  by  X/n^,  and  Xg  is  the  guide  wavelength. 

Typical  curves  of  the  frequeneyVversus  1/Xg  for  mode 
propagation  in  a  dielectric  waveguide.  Each  mode  in 
represented  by  a  line  which  is  confined  to  the  region 
between  the  lines  whose  slopes  are  c/n2  and  c/n^.  At 
the  frequency  1/  '  the  TEqi  mode  has  a  guide  wavelength 
of  Xg  ,  phase  velocity  Vp^  and  group  velocity  v^group. 

Field  plot  in  the  core  for  the  TEQ2  mode  far  from  cut-off 
and  for  a  small  difference  in  indices  of  refraction  of 
the  core  and  cladding. 

Field  plot  in  the  core  for  the  mo<*e  far  from  cut-off 

and  for  a  small  difference  in  indices  of  refraction  of  the 
core  and  cladding. 

Field  plot  in  the  core  for  the  EH.^  mode  far  from  cut-off 
and  for  a  small  difference  in  indices  of  refraction  of  the 
core  and  cladding. 

Field  plot  in  the  core  for  the  HE21  mode  far  from  cut-off 
and  for  a  small  difference  in  indices  of  refraction  of  the 


core  and  cladding. 


I .  INTRODUCTION 


In  a  light  pipe  electromagnetic  energy  is  propagated 
down  the  pipe  by  reflection  from  the  walls  of  the  structure. 

If  the  transverse  dimensions  are  comparable  to  the  wavelength  of 
the  light,  only  certain  field  distributions,  or  modes,  will 
satisfy  Maxwell's  equations  and  the  boundary  conditions.  In 
this  case  the  light  pipe  is  more  appropriately  considered  as 
a  waveguide.  Even  in  very  large  structures  there  are  waveguide 
modes,  but  there  are  so  many  of  them,  their  number  increasing 
as  the  area,  that  in  most  cases  a  geometrical  optics 
description  is  more  fruitfial. 

Waveguides  were  first  dealt  with  by  Lord  Rayleigh.^- 
Later  the  dielectric  waveguide  was  investigated  theoretically 
by  Hondros  and  Debye2  and  experimentally  by  Schriever.3 

The  distinction  between  metallic  ana  dielectric 
waveguides  is  in  the  reflection  mechanism  responsible  for 
confining  the  energy.  The  metallic  guide  does  so  by  reflection 
from  a  good  conductor  at  the  boundary.  In  the  dielectric 
waveguide,  this  is  accomplished  by  total  internal  reflection, 
which  is  gotten  by  having  the  central  dielectric  made  of  a 
material  of  higher  index  of  refraction  than  the  surrounding 
dielectric.  The  two  regions  will  henceforth  be  referred  to  as 


the  core  and  cladding. 
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In  a  metallic  guide  there  are  two  sets  of  solutions, 

the  transverse  electric  and  transverse  magnetic  modes.  In  the 

dielectric  guide  all  but  the  cylindricaily  symmetric  modes, 

TEQm  and  TM^,  are  hybrid,  i.e.  they  have  both  electric  and 

4,5 

magnetic  z-components .  In  general,  one  would  expect  two 
sets  of  such  hybrid  modes,  because  the  boundary  condi tions 
give  a  characteristic  equation  which  is  quadratic  in  the 
Bessel  functions  describing  the  field  in  the  central  dielectric. 
Beam  et  al6  gave  the  two  sets  for  n=i,  and  Abele7  arrived  at 
the  two  sets  by  a  graphical  solution  of  the  characteristic 
equation. 

Until  recently  the  main  concern  has  been  with  the 
three  lowest  order  modes,  HE^,  TEq^  and  TMq^.  However,  with 
the  increased  interest  in  end-fire  antennas  and  the  observations 
of  waveguide  modes  in  the  visible  region  of  the  spectrum,®  the 
higher  order  modes  assume  greater  importance. 

In  this  paper  th:  two  sets  of  solutions  are  given, 
including  cut-offs,  field  distributions,  and  conditions  far 
from  cut-off.  The  designations  of  the  hybrid  modes  are 
considered.  Some  properties  of  the  propagating  modes  are 


discussed. 
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II.  BASIC  EQUATIONS 

The  cylindrical  dielectric  waveguide  consists  of  a  core 
of  high  dielectric  constant  and  radius  a.  surrounded  by  a 
cladding  of  lower  dielectric  constant  Both  regions  are 

assumed  to  be  perfect  insulators  with  the  free  space  magnetic 
permeability  Such  a  structure  can  have  an  infinite  number 
of  modes,  but  for  given  values  of  6^,  and  £.*  only  a 

finite  number  of  these  are  waveguide  modes  have  their  fields 
localized  in  the  vicinity  of  the  core.  The  other  "unbound* 
modes  would  correspond  for  example  to  light  striking  the  core 
from  the  side,  passing  on  through  the  core  and  emerging  froui 
the  other  side. 

Choose  a  cylindrical  coordinate  system  r,  9,  z  with 
the  z-axis  lying  along  the  guide  axis.  A  waveguide  mode  is  a 
coherent  distribution  of  light,  localized  in  the  vicinity  of 
the  core  by  total  internal  reflection,  and  which  propagates 
down  the  guide  with  a  well  defined  phase  velocity.  That  is- 
the  z  and  time  dependences  are  given  by  exp  f i(hz-«t)}  , 
vdere  w  is  the  angular  frequency  and  h  is  the  propagation  constant 
which  is  determined  from  the  boundary  conditions. 

Because  of  the  cylindrical  symmetry,  the  other  components 
of  the  fiair'1  can  be  expressed  in  terms  of  Ez  and  Hz.9  The 
z-coroponents  of  the  field  satisfy  the  wave  equation  in  cylindrical 
coordinates.  The  solutions  are: 


r 


Ez  =  An  Jn  (?v1  r)  cos  (n9  +  f  n)exp  (i(hz-Wt)}  , 
H2  =  Bn  Jn  (A1  r)  cos  (n9  +  ^,n)exp  |i(hz-0fc:)}  . 


The  field  in  the  cladding  is  given  by  replacing  the  constants 
An  and  Bn  by  Cn  and  Dn#  and  by  replacing  the  Bessel  function 
Jn  (^l  r)  by  the  modified  Hankel  function  of  the  first  kind 
Kjj  (7v2  r)  •  These  particular  Hankel  functions  are  required  in 
the  cladding,  because  they  are  the  only  cylindrical  functions  that 
vanish  sufficiently  rapidly  as  r  increases  to  infinity  to  describe 
a  field  bound  to  the  central  dielectric.  With  the  definition 
of  the  propagation  constant  k2  =  w2|i£  ,  the  Tv's  are  defined  by 

\2  =  kx2  -  h2,  Tv22  »  h2  -  k22; 

both  are  real,  (p  n  and  y,  n  are  phase  factors,  which  are  related 
by  the  boundary  conditions. 

The  transverse  components  of  the  field  can  be  expressed 

,  g 

m  terms  of  E  and  H_  by 
z*  z 


E_  =  i 


p +  '&>  I  jL5z_1 
*-h2  l  b  x  hr  a  9  J 


Ee  -  i 


Hr  *  i 


He  -  i 


k2-h2 


_  [l  £2  _  b  Hz  1 

Lr  d8  h  d  r  J 

r  *2 1  jgz  +  i 

L .  n<tfh  r  cJ  9  b  r  J 


h  f  k2  b  E  1 

_  _  _ _ z.  +  _ 

k2-h2  Luu/h  b  r  r 


a  hz 
b  9  . 
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At  the  boundary  r  =  a,  the  continuity  of  the  tangential 
components  of  the  field  give  the  following  four  equations  for 
the  constants  An#  Bn,  C^,  Dn.  For  simplicity  7\-j_a  an^  ^2a  ^ave 
been  replaced  by  u  and  w.  The  continuity  of  the  tangential  components 
of  E  give 


Sii 

(iWu 


An 

J  = 
n 

=  C 

n 

V 

(8) 

A,. 

*1 

nh 

v* 

J 

n 

sin 

(n©  +  (prv) 

+  sn 

uU 

u 

'  cos  (n9  n) 

(9) 

= 

“Cn 

“T 

Wz 

I^sin  (n©  +  pn)~D:i 

w 

F^cos  (n©  +p  n) , 

the 

tangential  components 

of  H 

t 

Bn 

Jn  : 

=  Dn 

V 

(10) 

cos 

(n©  +p  n) 

-  B 

Bn  u2 

Jn 

sin 

(n©  ~^n) 

n  p.Ww 


Kn1  cos  (n©  +  <pn)  +  Dn  ^  kr  sin  (n©  +^n)  (11) 


The  primes  on  Jn  and  refer  to  differentiation  with  respect 
to  thpir  arguments,  u  and  w,  respectively.  Further,  define 
ni  and  n2  by 


nl  = 


u  Jn{u) 


n2 


w  (w) 


(12) 


For  the  four  equations  to  be  consistent,  the  determinant  of  the 
coefficients  must  vanish,  giving 

(n]_  +  n2)  Oc^  n^  +  k2  n2)  -  sin(n©  +p  n)  sin(n©  +^n) 

n2h2  (1/u2  +  l/w2)2  cos(n©  +p  n)  cos(n©  +^n) 


(13) 
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The  left  side  of  (13)  is  independent  of  angle?  therefore,  the 
phases  p  n  and  n  must  be  related  such  as  to  make  the  right 
s’de  a  constant.  This  will  be  the  case  if 

?n  ~Pn  =  (14) 

Then  (13)  becomes 

(n^  +  1*2)  (^1^  nl  +^22  n2^  =  n^h2(l/u2  +  1/w2)2.  (15) 

Eq.  (15)  together  with  Eqs.  (3)  determine  the  value  of  h. 

The  quantity  u  enters  (7)  both  explicitly  and  as  the  argument 
of  n-^.  However,  n-^  is  a  rapidly  varying  oscillatory  function  of  u. 
Hence,  (15)  can  be  considered  roughly  as  a  quadratic  equation 
in  n-^.  The  two  sets  of  solutions  are  the  two  sets  of  hybrid 
modes . 

From  the  set  of  four  boundary  condition  equations, 
the  relative  amounts  of  Ez  and  Hz  in  a  hybrid  mode  can  be 
found  in  a  straightforward  way.  The  result  is 

p  -  IW  Bn  cos (nO  +  ffn)  =  n(l/u2  +  1/w2)  (16) 

h  An  sin(n9  +  ^n)  nl  +  n2 

The  coefficient  \iuif. h  in  the  definition  of  ?  has  been  included 
for  convenience  in  later  use. 

Expressed  in  terms  of  P,  the  field  in  the  core  becomes 
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sz  =  Jn  (*1*>  FC' 

S,  -  i  -k-  I  Jn1  -  p  I  F 

r  *1  L  n  X^rJ  c/ 


Alr  J 

~  1  “x7  PJnL  ~  Fs 
1  L  XXr  J  ' 


_  _  ~h—  pj 

F 

n 

s 

Tc  2 

=  -i  JL L— 

P  h2  j  1  -  DH.n 

M^Al 

l  ki2  n  Xjr 

k  ^ 

=  i  ki 

Jnx  -  p  2£a 

pAA]_ 

"  *12  Xjr 

The  prime  indicates  differentiation  with  respect  to  the 
argument  of  the  Bessel  function  and 

Fc  =  An  cos  (n0  +^n)  exp  |  i(hz-att)J 

Fs  =  sin  (n0  +^n)  exp  {i(hz-cJt )j. 
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III.  THE  CUT-OFF  CONDITIONS 


The  cut-offs  for  the  various  modes  are  found  by 


s 

solving  (15)  in  the  limit  of  w*->  0. 


Define  new  quantities  <£  ^  and  2  by 


<1- 


J~-l 


By  using  the  Bessel  and  Hankel  function  identities,  Eqs.  (A4)  and 
(A6)  in  the  appendix,  (15)  can  be  rewritten  as 


Consider  first  the  case  of  n>2.  From  Eq.  (A10)  as 


w2-^  0, 


4  2  m  ^^(n-l)]  . 


The  two  sets  of  roots  of  (20)  are 

£  _  _I_  k?2 _ 

^  1  ~  n-1  kx2  +  k22  7 


-i  •>  “i 


^  is  a  function  of  u,  but  as  pointed  out  earlier,  varies 
so  much  more  rapidly  than  u  that  where  u  appears  explicitly 
in  (20)  it  can  be  considered  a  constant.  The  two  sets  of 
cut-off  values  for  different  n  can  be  obtained  by  substituting 
for  the  limit  of  w*-+0. 
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^  _  ki2  +  k22  i  ^ 

”  -  i^r- 

Eq.  (23)  is  equivalent  to  Jn(u)=0. 

The  cut-off  values  from  (22)  are  given  by 
Schelkunoff Those  from  (23)  were  found  by  Abele^  by  a 
graphical  solution  of  equation  (15)  .  Substituting  for 
in  terms  of  the  Bessel  functions  and  using  (Al)  with  n 
replaced  by  (n-1)  puts  Eq.  (22)  in  a  more  convenienc  form  for 
small  differences  in  dielectric  constant  between  core  and 
cladding.  The  result  is 


uJn-2<u>  ,  ,,  €l  "€.2 

- -  =  -  (n-1)  -  • 

Jn-1  ^  6  2 


(24) 


For  small  differences  in  dielectric  constant  the  cut-offs  are  given 

approximately  by  Jn_2 M  =  0.  For  n=2,  this  gives  another  set 

of  modes  whose  cut-offs  are  close  to  the  TE  and  TM  modes. 

om  om 

The  interesting  effects  this  has  on  optical  waveguide  modes  are 
discussed  in  a  companion  article  in  this  journal. 

For  n=l,  from  Eq.  (A9)  in  the  limit  of  w2-^0  the 
two  roots  of  (20)  are 


.  =  kl2  +  k22 

1 

1 

w2 

2  ko2 

1  k^2  +  k22 

In 

00 


(25) 


2 

i  w 


oo 


(26) 
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Both  of  the  above  equations  have  cut-offs  given  by 

Jx  (u)=0  (27) 

The  root  of  Eq .  (27)  at  u=0  corresponds  to  the  well-known 
HEn  mode  which  does  not  have  a  cut-off. 

Eq.  (27)  specifies  two  sets  of  modes  whose  cut-off 
conditions  are  identical.  However,  the  HE.^  mode  is  the  first 
mode  of  only  one  of  the  sets.  This  is  because  1/u2  in  equation 
(20)  cannot  be  considered  as  slowly  varying  near  u=0.  In  fact, 

£ ^  becomes  2/u2  near  u»0  and  the  quadratic  term  ,  drops 
out  of  the  equation. 

The  relative  amounts  of  Ez  and  Hz  in  the  hybrid 
modes  can  be  found  by  substituting  in  Eq.  (16)  the  limiting 
values  of  £  ^  and  £2  as  w-*»0.  The  results  together  with  the 
cut-off  conditions  are  summarized  in  Table  I. 

Let  u^  be  the  value  u  assumes  at  cut-off  for  the 
m'th  root  of  the  cut-off.  condition  involving  the  n'th  order 
Bessel  function.  The  possible  values  of  u^  are  the  roots  of 
the  equations  giving  the  cut-off  conditions  in  Table  I.  At 
cut-off,  w»0  and  h-k2;  hence,  the  first  of  Eqs.  (3)  becomes 

unm  =  2  If  (a A)  (ni2  -  n22)1/2  ,  (28) 

where  X  is  the  free  space  wavelength  and  n^  and  n2  are  the 
indices  of  refraction  of  core  and  cladding,  respectively. 

The  modes  which  can  propagate  are  those  for  which 

1/2 

u nrn  are  less  than  2  <a/x)  <nl2  ■  n22)  •  Since  u^  forms 
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an  increasing  sequence  for  fixed  n  and  increasing  m  or  for 
fixed  m  and  increasing  n,  the  number  of  allowed  modes  increases 
as  the  square  of  the  radius  a, 

IV.  CONDITIONS  FAR  FROM  CUT-OFF 

The  field  in  a  mode  far  from  cut-off  can  be  found  in 
the  same  way  as  that  near  cut-off.  Only  now  the  asymptotic 
forms  for  large  w  are  substituted  for  £2  in  E<3*  (20).  Table 
II  summarizes  the  conditions  far  from  cut-off. 

V.  MODE  DESIGNATIONS 

For  cylindrical  metallic  waveguides  the  modes  are 
designated  TE^  and  TM^,  or  in  the  older  literature  these 
were  and  E^  modes,  respectively.  The  transverse  electric 
mode  TE^  can  be  derived  from  a  single  field  quantity,  the 
z-component  of  the  magnetic  field,  hence  the  alternative 
designation  for  this  mode. 

In  the  dielectric  waveguide  only  the  cylindrically 
symmetric  n  =  0  modes  are  either  transverse  electric  (TEojn) 
or  transverse  magnetic  (™om)  •  The  otner  modes  are  hybrid, 
i.e.  they  have  non-zero  values  for  both  Ez  and  Hz . 

Following  Wegner,  10  Beam^  suggested  a  scheme  for  the 
designation  of  the  hybrid  modes  based  on  the  relative 
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contributions  of  E2  and  Hz  to  a  transverse  component  of  the  field 
at  some  reference  point.  If  Ez  makes  the  larger  contribution, 
the  mode  is  considered  E-like  and  designated  EH^,  etc.  This 
method  of  designation  is  arbitrary,  for  it  does  depend  on  the 
particular  transverse  component  of  the  field  chosen, the  reference 
point,  and  how  far  the  wavelength  is  from  cut-off.  However, 
the  use  of  two  letters,  such  as  EH  and  HE,  to  designate  the 
hybrid  modes  is  reasonable  because  it  does  imply  the  hybrid 
nature  of  these  modes. 

It  has  become  common  usage  in  the  microwave  literature 
to  refer  to  the  mode  without  a  cut-off  as  the  HEj^  mode. 

Referring  to  Table  II,  it  is  seen  that  this  mode  has  a  value 
of  P  =  -1  far  from  cut-off.  It  is  proposed  that  all  the  modes 
with  P  =  -1  be  designated  HE^  and  the  modes  with  P  =  +1  be 
designated  EH^.  The  basic  physical  difference  between  the  HE 
and  EH  modes  will  be  discussed  later  in  the  section  on  field 
plots . 

The  subscripts  on  KE^  or  EH^  refer  to  the  n ’  th 
order  and  m*th  rank,  where  the  rank  gives  the  successive  solutions 
of  the  boundary  condition  equation  involving  Jn.  It  is 
customary  to  label  these  solutions  in  order  starting  from  m  •  1; 
this  procedure  is  followed  here.  Notice  that  this  makes  the  cut-off 
parameters  for  the  HEj^  and  EH^  mode  the  same,  namely  th».  ~t  of 
J^tu)=0  at  u=3.832.  Table  I  summarizes  the  mode  designations  su..^es.'' 


here . 
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VI.  THE  TRANSVERSE  COMPONENTS  OF  THE  FIELD 


By  use  of  Eqs.  (Al) ,  (A2)  and  (17)  the  field  in  the 
core  can  be  written  as 


E„  = 


= 


En  - 


Hz  -  - 


K  = 


2  Jn-1  2  Jn+1  J  Fc' 


Jn  <*!*>  Fc' 

.  _h_  [  Izl 

1  Xi  L  2 

P  Jn  Fs* 

-i  ki2  L  1  ~ 


n-i  2 


1  * 

»* r  —  i 


F 


s ' 


h 


i_  J. 


(29) 


n-1 


1  +  Ph^/kj2 


He  *  i 


-  PhVkj2  j  _  1  +  ?h2/k, 2 


'n-1 


J 


n+1 


l]Fs- 

]'«* 


The  quantities  Fc  and  Fg  continue  to  be  given  by  Eqs.  (18). 

For  small  differences  in  refractive  index  between  the  core 

and  cladding  ki2/k2^  =  1,  h  /k^2  =*  1  and  from  Tables  I  and  II 

P  *  +  1.  The  value  P  =  +1  is  for  the  EH  modes  and  P  =  -1  goes 

with  the  HE  modes.  It  is  seen  from  Eqs.  (29)  that  the  transvers 

components  of  the  field  depend  on  r  through  J  for  the  EH 

n+1 

modes  and  have  a  dependence  of  for  the  HE  nodes. 


» 
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VII .  ENERGY  FLOW 


The  energy  flow  per  unit  area  as  a  function  of  r  and  9 
is  given  by  the  real  part  of  the  complex  Poyntir.g  vector, 

S  =  1/2  E  x  H,  where  H  is  the  complex  conjugate  of  H. 

4M  ««*  AM  mm 

Only  the  z-component  of  J3*  is  real  and  the  power  per  unit  area 
is  given  by 

Sz  =  1/2  (Er  He  -  Ee  Hp)  .  (30) 

Substitution  from  Eqs.  (29)  gives  af*-er  some  simplification 

2 


i  h_  ki?  »  2  f  (1-P)  (l-Ph2AV 
=  2  An  [  4 


) 


'n-l 


+  (l+p)  (l+*h2Ai2)  jn+12 


(31) 


2v2  a,_  2 


_  X"P  h  A1 


Jn-1  Jntlcos2^ne+^ 


'*e+^n)j 


For  small  index  differences  P  =  +  1  and  the  term  containing 
0  drops  out.  The  energy  flow  is  then  nearly  circularly 
symmetrical  with  a  radial  dependence  of 

sz  *  Jn5a  2  (32> 

The  upper  sign  n-l  is  for  the  HE  modes  and  the  lower  for  EH. 


i 


vWxv.AU 
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VIII.  CHARACTERISTICS  OF  FROPAGATIN'G  MODES 

For  a  waveguide  mode  and  must  both  be  real. 

Hence  from  Eq.  (3) 

k22  <  h2  <  k]2.  (33) 

A  guide  wavelength  kg,  phase  velocity  vp^  and  effective 
refractive  index  can  be  defined  for  a  mode  at  the  free 

space  wavelength  A  by 

W/n  =  vph  =  =  c/r.eff,  (34) 

whe’  e  c  is  the  velocity  of  light  in  vacuum..  With  (34)  and  the 
definition  of  the  k’s,  Eq.  (33)  becomes 

c/nj_  £  vp’n  <  c/n2  »35) 

That  is,  the  phase  ’velocity  is  intermediate  between  the 
velocities  of  propagation  in  bulk  mate  'i?.l  of  which  the  core  ana 
the  cladding  are  made.  At  cut-off  h»k2  or 

vph  =  c/n2*  (36) 

Far  from  cut-off  the  ether  equality  in  (35)  holds. 

Eq.  (36)  is  the  physical  optics  analogue  of  the 
geometrical  optics  condition  that  propagation  in  the  light  pipe 
occurs  when  the  angle  of  incidence  or.  the  boundary  exceeds  the 
critical  angle  for  total  internal  reflection.  Consider  a  plane 
wave  incident  on  the  side  of  the  core  as  shown  in  figure  1. 

The  wave  normal  is  in  the  direction  p,  and  S.  and  S2  are  two 
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equiphase  surfaces  separated  by  X/n^,  the  wavelength  in  a 
homogeneous  medium  of  refractive  index  nn .  The  apparent 
wavelength  along  the  guide,  Xg,  is  given  by 

Xg  =  Vni  sin  .  (37) 

But  Snell's  law  gives  for  the  critical  angle  for  total  internal 
reflection,  n2  =  sin  <X  -  Hence  at  cut-off  (37)  becomes 
Xg  =  X/n2,  or  Vpk  =  c/n2,  which  is  the  same  as  Eq.  (36). 

Fig.  2  gives  the  very  useful  plot  of  frequency  versus  the 
inverse  of  the  guide  wavelength,  that  is  c/X  vs.  1/Xg  for  the 
propagating  modes.  To  obtain  these  curves,  Eqs.  (3)  and  (15) 
are  solved  for  h  as  a  function  of  X  and  the  parameters  of  the 
guide  for  each  propagating  mode.  The  mode  lines  shown  in  Fig.  2 
are  only  schematic.  A  number  of  machine  computations  of  these 

lines  are  available  in  the  literature  for  the  HE-^  and  TEqi 

,  6,11 

modes , 

Propagation  in  the  modes  is  such  as  to  be  confined  to  the 
region  between  the  lines  whose  slopes  are  c/n^  and  c/n2 . 

Each  mode  as  a  function  of  wavelength  is  represented  by  a  line  which 
approaches  the  c/n^  line  far  from  cut-off  and  terminates  at  the 
c/n2  line  at  cut-off.  All  the  modes  which  have  cut-offs 
terminate  sharply  at  the  c/n2  line,  but  the  HE^  mode,  which 
does  not  have  a  cut-r.f,  approaches  the  c/n2  line  slowly,  finally 
merging  with  it  at  the  origin.  The  number  of  modes  increases 
as  the  square  of  the  frequency. 
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The  significance  of  a  mode  line  can  best  be  understood 
by  a  specific  example.  Consider  the  TEq^  mode  excited  at  the 
frequency  V  ' .  Then  1/Ag  is  the  coordinate  of  the  intersection 
of  the  TEq^  line  with  the  ordinate  V'  .  The  slope  of  the  line 
connecting  this  point  with  the  origin  is  c/n^  (A) ,  where  n0^(A) 
is  the  effective  index  of  refraction  of  the  TEQ1  mode  excited 
at  the  free  space  wavelength  A.  c/nr>^(A)  is  the  phase  velocity, 
whereas  the  slope  of  the  mode  curve  at  V'  is  the  group  velocity. 
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IX.  FIELD  PLOTS 

The  field  distribution  can  be  given  in  the  usual  way 
by  field  lines  in  which  the  direction  of  the  line  at  a  point 
gives  the  direction  of  the  field  and  the  density  of  lines  its 
magnitude.  Only  the  field  in  the  core  will  be  considered.  The 
field  in  the  cladding  can  be  inferred  from  the  boundary  conditions 
which  require  that  the  tangential  components  of  E  and  H,  and  the 

AS*  M* 

normal  components  of  £E  and  (iH  be  continuous.  Since  the  field 

A/^  AM 

components  in  the  cladding  depend  on  modified  Hankel  function 
which  monotonously  go  to  zero  with  increasing  r,  the  density 
of  field  smoothly  decreases  with  an  increase  in  the  radial 
coordinate . 

In  Fig.  3  is  shown  the  field  distribution  in  the  core 
for  the  TE^  mode  far  from  cut-off.  The  field  components  from 
which  Fig.  3  was  sketched  are  given  in  Eqs.  (17)  or  (29) . 

The  distribution  is  the  instantaneous  values  in  a  transverse 
plane  and  in  two  longitudinal  half  planes  of  length  Xg/2  in 
the  z-direction.  In  the  transverse  plane  the  magnetic  field  lines 
shown  do  not  form  closed  curves;  this  is  done  to  imply  that  these 
field  lines  close  by  going  down  the  z-direction.  The  dots  and 
crosses  in  the  longitudinal  planes  specified  that  the  electric 
field  is  perpendicular  to  these  planes.  The  field  line  enters 
the  plane  at  a  dot  and  leaves  it  at  a  cross. 


-19- 


The  z-component  of  the  magnetic  field  varies  as 
Jo^u  r/a) .  From  Table  II  far  from  cut-off  u  is  the  second  zero 
of  J^u)  =  0  or  u  =  7.02,  and  at  cut-off  u  is  the  second  zero 
of  Jq  (u)  =  0  or  u  =  5.52.  Hence,  as  the  cut-off  wavelength 
is  approached  the  field  distribution  in  the  core  readjusts  itself 
so  as  to  make  the  boundary  of  the  core  shift  from  the  value 
u  =  7.02  shown  in  Fig.  3  to  the  value  u  =  5.52. 

For  the  TE0-|_  mode  the  boundary  r  =  a  is  at  (j.  =  3.83 
far  from  cut-off  and  at  |i  =  2.41  at  cut-off.  The  field 
distribution  for  the  TMqj  mode  is  obtained  from  Fig.  3  simply 
by  interchanging  the  roles  of  E  and  H. 

Since  the  Poynting  vector  is  a  function  of  the 
vector  product  of  the  transverse  components  of  the  field,  it  is 
clear  from  Fig.  3  that  an  image  of  the  intensity  distribution 
in  the  TEq2  or  TMq2  modes  should  consist  cf  two  concentric 
circles . 

It  is  interesting  to  compare  the  field  distribution 

for  the  TEQm  and  TMQin  modes  with  that  obtained  for  the  metallic 

12 

waveguide.  In  the  later  case  the  electric  field  is  normal 
to  the  metallic  boundary  surrounding  the  core  and  the  magnetic 
field  is  parallel  to  it.  Hence,  the  metallic  TEom  modes 
look  like  the  dielectric  TEoro  far  from  cut-off,  but  the  metallic 
TMojn  modes  have  the  field  distrubition  cf  the  dielectric 
Tfl  modes  at  cut-off. 


9 


For  the  hybrid  modes  (n  ^  1}  the  field  distribution 
simplifies  considerably  for  the  case  of  a  small  index  difference 
between  the  core  and  the  cladding.  From  Eq .  (29)  for  the  z-components 


given  by 


(Ez,  H  )  (X  j  (ur/a)  cos  n9, 
"  "  n 


the  transverse  components  of  the  electric  field  are 
Er  c<  (±1)  (u  r/a)  cos  n9, 


Eg  ©<  Jn+i  (u  r/a)  sin  n0. 


The  proportionality  constants  are  the  same  for  Er  and  E^.  The 

plus  sign  is  for  the  EH^  modes  and  the  minus  sign  for  the 

HE  modes, 
nra 

A  given  field  line  will  be  contained  on  a  surface 
whose  projection  in  the  transverse  r,  9  plane  is  the  solution  of 


=  Er/Ee 


For  a  small  difference  in  refractive  indices  between  the  core 
and  the  cladding  Eqs.  (39)  and  (40)  apply,  and  the  right  side 
of  (41)  is  a  function  of  0  only.  The  result  is 


dr 
r  d0 


+  cos  n  9 
sin  n  9. 


The  above  can  be  integrated  to  give 
r  ®  C  (sin  n  9)  +.  1/n, 
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where  C  is  an  integration  constant.  The  plus  sign  in  the 
exponent  is  for  the  EH^  modes  and  the  minus  sign  for  the 
HEnjn  modes.  By  assigning  different  values  to  C  a  family  of 
curves  is  obtained  which  gives  the  electric  field  in  a  transverse 
plane.  The  magnetic  field  can  be  found  in  the  same  way; 
it  is  the  same  as  for  the  electric  field  but  with  the  pattern 
rotated  by  T?/(2n)  . 

For  the  HEj^,  modes  Eq .  (43)  becomes 
r  sin  0  =  constant. 

This  is  the  equation  for  a  set  of  straight  lines  parallel  to  the 
x-axis.  Fig.  4  shows  the  field  distribution  for  the  HE-j^ 
mode.  Far  from  cut-off  the  boundary  of  the  core  is  at  the  second 
root  of  JQ(u)  =  0  or  u  =  5.52  and  at  cut-off  the  boundary  is  at 
the  zero  of  J^(u)  =  0  at  u  =  3.83.  For  the  HE-}  ^  mode  the 
boundary  shifts  from  u  «  2.41  for  short  wavelengths  to  u  =  0 
with  increasing  wavelength. 

From  Eq.  (43)  the  electric  field  lines  in  a  transverse 
plane  for  the  EH^  mode  satisfy 

r  =  C  sin  9,  (45) 


or 

x2  +  (y  -  C/2)2  =  C2/4. 


(46) 


By  assigning  different  values  to  C,  Eg.  (46)  gives  a  set  of 
circles  that  are  all  tangent  to  the  x-axis  at  che  origin.  Fig.  5 
gives  the  field  plot  for  the  EH^  mode.  Far  from  cut-off  the 
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boundary  of  the  core  is  at  the  first  zero  of  (u)  =  0 

or  u  =  5.14  and  at  cut-off  the  boundary  is  at  the  first  zero 

of  J^(u)  =  0  at  u  =  3.83  . 

For  the  HE21  mode  the  electric  field  in  a  transverse 

plane  is 

r  =  C  (sin  20) ~1/2  (47) 


or 


xy  =  Constant.  (48) 

For  various  values  of  the  constant  2q.  (48)  gives  a  set  of 
hyperbolas.  The  field  lines  are  shown  in  Fig.  6.  Far  from 
cut-off  the  boundary  of  the  core  is  at  the  first  root  of 
Jx(u)  =  0  or  u  =  3.83  and  at  cut-off  the  boundary  is  at  a  value  of 
u  somewhat  larger  than  the  first  root  of  JQ(u)  =  0  or  u  =  2.41. 

To  find  a  more  precise  value  of  u  at  cut-off  for  the 
case  of  a  small  index  difference  the  left  side  of  Eq.  (24)  can  be 
expanded  about  the  zero  of  Jn-2 =  0.  To  the  first  order 
terms  the  result  is 


u ' 
nm 


u 


n-2  ,m 


n-1 


un-2. 


m 


(49) 


The  prime  on  u^  indicates  that  it  is  the  cut-off  parameter 
for  the  HE^  modes  obtained  from  Eq .  (24). 

Where  the  field  lines  in  Fig.  6  do  not  form  closed 
curves  but  end  in  a  plane,  the  field  lines  close  by  moving 
perpendicular  to  the  indicated  plane.  The  lower  portion  of 
Fig.  6  shows  the  closed  contours  for  the  electric  field  in  a 


hyperbolic  section  bb' . 


ft  ■■  n'fww;, »  f  "MO"  vj»mr»r.  f>  fcJg^,'^5^ivU4JW>i!f  uamgqwsnf 


TT7*TT 
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The  field  plots  for  the  higher  order  modes  are 

obtained  in  the  same  way  as  those  above. 

The  EH^  and  HE^  modes  both  have  roughly  the  same 

r-dependences  for  (Ez,  Hz) ,  mainly  Jn(ur/a).  *But  the  transverse 

components  depend  on  Jn+^(ur/a)  for  the  EH  modes  and  on 

Jn_^(ur/a)  for  the  HE  modes.  This  means  that  the  field  lines 

which  are  parallel  to  the  guide  axis  at  z=0  tend  to  form  closed 

contours  by  going  to  larger  radii  in  the  region  0  <|z|<Xg/4 

for  the  EH  modes,  but  on  the  whole  close  by  going  to  smaller 

radii  for  the  HE  modes.  Hence,  for  the  EH  modes  the  peaks  in 

the  Poynting  fluxare  located  further  from  the  center  of  the 

guide  than  the  peaks  in  E  and  H  •  the  reverse  is  the  case  for 

z  z 

the  HE  modes. 
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Table  I.  Summary  of  Cut-Off  Conditions 

The  Bessel  function  of  order  n  and  argument  u  is 
given  by  Jn(u),  and  ^  and  6  ^  are  t^ie  dielectric  constants 
for  core  and  cladding.  P  gives  the  relative  amount  of 
H  to  E-  in  a  mode  (see  text  for  exact  definition) . 
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Table  II.  Summary  of  Conditions  Far  from  Cut-Off. 


First  Solutions 

Second  Solutions 

u 

P 

u 

P 

n**0 
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Jn-l(u>=° 

-i  (HE) 

Jn+i(u)=0 
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APPENDIX:  PROPERTIES  OF  BESSEL  AND  HANKEL  FUNCTIONS 

The  argument  of  the  Bessel  functions  is  u  and  of  the 
modified  Hankel  functions  w. 


n  Jn  =  1  „ 

2  (Jn-l  +  un+l 


u 


). 


(Al) 


V 


-n 


*  i  {Jn-l  "  Jn+1> ' 
2 

-  (-l)n  JR. 


(A2) 


From  the  first  two  equations  above 


n 

u  Jn  u  Jj 


*  °n-l  -  n 


(A4) 


u1 


For  w  real  the  modified  Hankel  functions  K  (w)  are  defined  by 

n 


Kn(w)  =  JL  i  n+1  Hn  (1>  (iw), 


(A5) 


where  H  ^  (iw)  are  the  Hankel  functions  of  the  first  kind.  The 
n 

equation  corresponding  to  (A4)  is 


V  »  *n-l  + 

w  Kn  w  Kn  w2 


(A6) 
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For  small  w. 


K0(w) 

=  In  (2/^T  w) , 

(A7) 

V"> 

=  (n-1);  2n”^  n~2  for  n>l. 

(A8) 

where  ?  is  Euler ' 

s  constant  and  equal  to  1.781. 

Still  in  the 

limit  of  w  small, 

= 

W  Kx 

In  (2/2Tw), 

(A9) 

V-l  = 

w  K„ 

j^2  (n-1)  J  1  for  n— 2. 

(A10) 

The  asymptotic  expressions  for  large  w  are 

K  (w) 
n 

=■  ( 7 T  /2w)  exp  \  -w  |  , 

(All) 

*n-l  _ 

1 

(A12) 

wKn 

w 

CAPTIONS  FOR  FIGURES 


Construction  to  show  the  equivalence  of  the  waveguide 
condition  that  at  cut-off  v  ^  =  c/n2  with  the 
geometrical  optics  conditions  that  propagation  occurs 
only  if  the  angle  of  incidence  of  the  wave  on  the 
fiber  wall  exceeds  the  critical  angle  for  total 
internal  reflection.  The  wave  normal  is  given  by  p, 

S^  and  S2  are  two  equiphase  surfaces  separated  by 
A/n^,  and  Ag  is  the  guide  wavelength. 

Typical  curves  of  the  frequency i'  versus  1/Ag  for  mode 
propagation  in  a  dielectric  waveguide.  Each  mode  is 
represented  by  a  line  which  is  confined  to  the  region 
between  the  lines  whose  slopes  are  c/n2  and  c/n^. 

At  the  frequency  V  '  the  TEq^  mode  has  a  guide 

wavelength  of  Ag',  phase  velocity  v  ^  and  group 
velocity  vgroup. 

Field  plot  in  the  core  for  the  TEg2  mode  far  from  cut¬ 
off  and  for  a  small  difference  in  indices  of 
refraction  of  the  core  and  cladding . 

Field  plot  in  the  core  for  the  HE12  mode  far  from  cut¬ 
off  and  for  a  small  difference  in  indices  of  refraction 
of  the  core  and  cladding. 
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l 

| 

I  Fig.  5  Field  plot  in  the  core  for  the  EHn  mode  far  from 

f 

j  cut-off  and  for  a  small  difference  in  indices  of 

! 

refraction  of  the  core  and  cladding. 

;  Fig.  6  Field  plot  in  the  core  for  the  HE21  mode  far  from 

t 

■'  cut-off  and  for  a  small  difference  in  indices  of 

5  "  refraction  of  the  core  and  cladding. 
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